Periodic Oscillations of Josephson- Vortex Flow Resistance in Bi 2 Sr 2 CaCu 2 
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To study the Josephson-vortex system in the intrinsic Josephson junctions, we have measured the 
flow resistance as a function of magnetic field parallel to the ab plane in Bi2Sr2CaCu2 08+i/ single 
crystals. Although it was known that the flow resistance increases smoothly with increasing field, 
we have found novel oscillations of vortex flow resistance in the small current regime. The period 
of the oscillations corresponds to the field which is needed to add "one" vortex quantum per "two" 
Josephson junctions. Commensurability between the lattice spacing of Josephson vortices along the 
ab plane and the size of the junction is related to the flow velocity of vortices. The results show 
that Josephson vortices form triangular lattice in the state where the oscillations occur. 

PACS numbers: 74.60. Ge, 74.25.Dw, 74.72.Hs, 74.25. Fy 



X 



When an external magnetic field is applied along the 
plane of Josephson junctions in the superconducting 
state, the field penetrates into the junctions as quan- 
tized Josephson vortices. Since high-T c cuprates such as 
Bi 2 Sr 2 CaCu208+; / (BSCCO) naturally consist of atomic- 
scale Josephson junctions (intrinsic Josephson junctions; 
IJJ's) based on the layered crystal structures and large 
anisotropy , the magnetic field penetrated along super- 
conducting planes forms a Josephson-vortex system. Re- 
cently, the Josephson-vortex system in high-T c cuprates 
has attracted much attention because of fundamental in- 
terests and possibilities for applications. 

When the intrinsic pinning by the superconducting 
layers [Q] becomes strong enough, the commensuration 
of the vortex lattice with the underlying layered struc- 
ture becomes important, so that the drastic changes in 
the vortex phase diagram is expected. To clarify ther- 
modynamic phases of the Josephson-vortex system, both 
theoretical and experimental studies have been done ex- 
tensively. Theoretically, several new phases [0, a smec- 
tic phase M, and a tricritical point where the order of 
the melting transition changes || have been predicted. 
Experimentally, the important facts such as the melt- 
ing transition of the Josephson-vortex systems H ||, 
the oscillatory melting temperature [9|, and the exis- 
tence of the vortex smectic phase M have been found 
in YBa 2 Cu 3 07-5 (YBCO). 

Interesting phenomena appear when the vortex dy- 
namics is influenced by a periodic one-dimensional pin- 
ning potential. Periodic oscillations of pinning strength 
as a function of the magnetic field H have been observed 
in an artificially structured superconducting film with 
channels of weak pinning in a strong-pinning environ- 
ment | JTc| . The commensurability between vortex lattice 
and the channel width is related to this phenomenon. 
Additionally, Gurevich et al. have reported the periodic 
oscillations of in-plane resistance R(H) in superconduct- 
ing multilayers of NbTi/Cu when the Josephson vortices 
move perpendicular to the plane [lll]. They concluded 



that the origin is a dynamic matching effect in terms 
of the generation of standing electromagnetic waves in 
moving vortex lattice along the c-axis. 

Although the influence of the periodic pinning po- 
tential is expected to be more important in highly 
anisotropic materials such as BSCCO than in YBCO and 
others, there is little information on the Josephson vor- 
tices in BSCCO. In this Letter, we report on novel peri- 
odic oscillations of the flow resistance of Josephson vor- 
tices as a function of the in-plane magnetic field. The 
period of the oscillations suggests that the ground state 
of Josephson vortices is in triangular lattice when the 
oscillations occur. 

High quality single crystals of BSCCO were grown by 
traveling-solvent floating-zone technique A platelet 
of single crystals was carefully cut into narrow strips. Af- 
ter four contacts with electrodes were formed by using a 
silver paste, the center of strips of the single crystals was 
cut by focused ion beam (FIB) . The same process to fab- 
ricate IJJ's has been used in the case of BSCCO whiskers 
fl3f . An FIB image of a typical sample is shown in Fig. 
1. The inset is the schematic illustration of the junc- 
tion. The dimensions (width (w)x length (I) x thickness 
(t)) of the measured samples A, B and C are 7.3x13x1.3, 
18x16x1.6 and 31x17x2.0 ^m 3 , respectively. The su- 
perconducting transition temperature T c is 86 K in all 
samples. 

While resistance of strips was smaller than lil before 
the fabrication of junctions, the resistance of sample A, 
B and C increased to 1400, 560 and 315 at room tem- 
perature after the fabrications, respectively. Therefore, 
the measured resistance is almost c-axis resistance of the 
junction parts and the contribution of in-plane resistivity 
can be ignored. In the small mesa crystals of BSCCO, 
the flow of Josephson vortices, which is not affected very 
much by pancake vortices, has been realized in the so- 
called lock-in state [p|-|l6|. The ab plane of the sample 
was possible to orient parallel to the magnetic field to 
an accuracy of 0.005 degrees by adjusting the flux-flow 
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resistance as a function of the misorientation angle to its 
maximum. 

When the field is applied parallel to the ab planes in 
the superconducting state, the current along the c-axis 
drives Josephson vortices to the direction perpendicular 
to both current and field. The motion of the vortices 
causes the flow resistance. Up to the present it has been 
reported that the flow resistance increases smoothly with 
increasing field because the number of vortices increases 
[ p5| , ^6[ . However, we have discovered novel oscillations as 
a function of field as shown in Fig. 2. When the magnetic 
field is increased, finite flow resistance appears in 2 ~ 6 
kOe depending on the sample and the current density. 
This field may be related to the pinning effect of Joseph- 
son vortices due to the quality of the surface boundary 
of samples. The oscillations start from a constant field 
-Hstart ~ 7-8kOc in all samples. The period H p of the 
oscillations is quite constant in the wide range of mag- 
netic fields. As the magnetic field is increased, the flow 
resistance increases accompanying with the oscillations 
and suddenly drops above a field H stop . 

ii/stop decreases with increasing the absolute value of 
the angle (8) between the magnetic field and the ab plane. 
This may be due to the breakdown of the lock-in state 
because the pancake vortices penetrate into the junc- 
tion and are pinned to stop the flow of Josephson vor- 
tices when the c-axis component of the magnetic field 
exceeds the c-axis lower critical field H c \. Actually, H st0 p 
is roughly proportional to 1/ sin 8 and the coefficient is 
of the order of H c \. The maximum of H st0 p could be 
reached to 40 kOe with our alignment set-up. On the 
other hand, .ffstart and H p does not almost change even 
if the field is intentionally shifted from the optimal posi- 
tion. 

The temperature dependence of flow resistance has also 
been measured from 4.2 K to T c . The periodic oscillations 
can be observed at all temperatures below (T c — 5) K. H p 
is independent of temperatures. These facts show that 
the oscillations relate only to the field component parallel 
to the ab plane, namely, the dynamics of the Josephson 
vortices. 

It is interesting that periodic oscillations can be ob- 
served only in sufficiently small current regime. I-V char- 
acteristics of the sample B are shown in Fig. 3. In the 
fields where the resistance is minimum, the nonlinear- 
ity of I-V curves becomes large and a kink structure is 
observed around 50/iA. The periodic appearance of this 
structure as a function of fields corresponds to the oscil- 
lations of flow resistance. In the larger currents > 200/iA 
the oscillatory behavior disappears thoroughly. 

H p depends on the dimensions of the junctions. To ex- 
tract H pi the power spectral densities for three samples 
which are obtained from the data in Fig. 2 by FFT are 
demonstrated in Fig. 4. The sharp fundamental peaks in- 
dicate that Hp is uniform over the wide fields. The small 
peaks corresponding to the second harmonic waves are 



observed because the oscillations are distorted slightly 
from sinusoidal. The positions of the peaks, which are 
inversely proportional to H p , are plotted as a function of 
the width (w) of the junctions in the inset of Fig. 4. We 
see that the l/-ff p is proportional to the width (w). 

What is the origin of these oscillations? It should 
be checked on the relation between H p and the num- 
ber of vortices in the junctions. Supposing that the 
Josephson vortices uniformly penetrate into every junc- 
tion, the increment of the field (Hq) which is needed 
to add one vortex per a junction is 4>q/ws, where 4>o 
(~ 2.07 x 10~ 7 G cm 2 ) is flux quantum and s (=15 A) the 
junction thickness, namely, the distance between super- 
conducting Cu02 layers. 1/Hq(w) is shown as a broken 
line in the inset of Fig. 4. Interestingly, the experimen- 
tal values of 1 /H p are close to the twice of 1/H (w) as 
shown by a solid line, which means that H p corresponds 
to adding "one" flux quantum per "two" IJJ's. We can 
check whether the resistance takes a maximum or a mini- 
mum when the field is nH p (n : integer) . Experimentally, 
we have shown that the oscillations have minima at nHp. 

The above experimental results can be explained by 
assuming that; (1) the lattice structure of the Josephson- 
vortex system is the triangular lattice, and (2) at both 
right and left surfaces of the junction, a surface barrier, 
i.e., Bean-Livingston barrier JT^ works on vortices, al- 
though both attractive image force and repulsive force 
by Meissner current are quite small because of a little 
inhomogeneity of the magnetic field around vortices in 
junctions. For example, let's consider that the field is 
n(Ho/2) nH p . The Josephson vortices are distributed 
as schematically shown in Fig. 5(a). This ordered state 
fits to the width of junctions. When the vortex lattice 
flow from left to right, vortices at right side exit beyond 
small surface potentials and vortices would enter at left 
side simultaneously. Now, the total potential which the 
lattice feels at both sides becomes maximum. At this 
situation, the averaged velocity and hence resistance of 
vortex flow has minimum. In a short time scale, the 
velocity of the vortex lattice would oscillate with time 
whenever vortices go over the boundaries. As the field 
is increased, additional vortices are forced to enter into 
the ordered lattice. The Josephson vortex lattice could 
move more easily because of the incommensurability be- 
tween the width of the junctions and the lattice spacing 
of vortices. Then, the flow resistance becomes larger than 
that in the above fitted case. When the field goes over 
(n + l/2)(i?o/2), the lattice starts to return to the fitted 
state. At (n+l)(Ho/2), matching between the Josephson 
vortex lattice and the size of junction recovers as shown in 
Fig. 5(b), where the Josephson vortices enter in alternate 
layers comparing with Fig. 5(a). Hence, it can be plau- 
sible that the oscillation period of the flow resistance is 
a half of Hq- When the Lorentz force from the applied 
current becomes large enough comparing with the sur- 
face potential, it is expected that the oscillations smear 
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out in this scenario, which can explain the I-V charac- 
teristics qualitatively. More complicated models are also 
possible, which may explain the details of flow dynamics. 

There are several theoretical studies with respect to 
the thermodynamic phase diagram of the Josephson vor- 
tex system in anisotropic layered superconductors 
Ivlev et al. Q have derived the structure of Josephson- 
vortex lattice as the compressed hexagons of triangular 
lattice pointing along the c-axis in high enough fields by 
using London theory. However, in the case of strong 
magnetic fields of order 4>o/js 2 , where 7 is anisotropy 
ratio, the three-dimensional(3D) anisotropic London or 
Ginzburg-Landau approach is not adequate This 
field is estimated as 46 kOe in 7 ~ 200 . It is delicate 
whether the London theory can be applied at the field 
where the oscillations occur. On the other hand, the ex- 
tensive Monte Carlo simulations with the 3D anisotropic 
frustrated XY model by Hu et al. have also shown that 
a triangular lattice is the ground state where the Joseph- 
son vortices distribute in every layer for the anisotropy 
parameter 7 >8 when the filling factor is 1/32 ||. For 
7 =150 as is the case of BSCCO, it is known that the 
compressed triangular lattice is achieved even below the 
tricritical point value (~ 15 kOe) by rescaling 0. These 
theoretical results are consistent with our experimental 
ones. For the magnetic fields lower than 7kOe, the os- 
cillations cannot be observed in all samples in spite of 
the existence of finite flow resistance. A reason why the 
oscillations smear out at the low magnetic fields may be 
that the distance between the vortices is not close enough 
to form a triangular lattice. Ivlev et al. have predicted 
that a triangular array becomes unstable with respect to 
the formation of several new types of vortex lattices in 
the low fields ||. In addition to the ground state, the 
existence of several metastable states has been indicated 
. In this situation, the periodic oscillations could not 
be expected. 

Similar phenomena have been reported in artificially 
structured materials as described above jl0|,[ll| , although 
the origins are different from our case. Recently, in 
YBCO single crystals of 60 K phase which is more 
anisotropic than 90 K phase, an oscillatory behavior of 
the in-plane resistivity in the field parallel to the ab plane 
have been reported |^|. In their experiments, the oscilla- 
tions have minima when the period of the layered struc- 
ture of YBCO matches the intervortex spacing in the di- 
rection normal to the layers. The period of the oscillation 
is proportional to B~ 5 and is independent of the size of 
samples. This phenomenon originates from the commen- 
surability between s and vortex-lattice space along the 
c-axis, while the two lengths, the junction width w and 
vortex-lattice spacing along the ab plane are important 
in present case. 

In summary, we have measured the flow resistance of 
Josephson vortices in IJJ's of BSCCO single crystals. 
The novel oscillations of the resistance as a function of 



field by the Josephson- vortex flow have been found. The 
period of the oscillations depends only on the width of 
junctions. The fact that the period is Ho/2, which is 
needed to add "one" Josephson vortex per "iW layers, 
suggests that the lattice structure of the Josephson vor- 
tices is triangular. The oscillations originate from com- 
mensuration between the triangular lattice and the width 
of the junction perpendicular to the field. At present, 
our measurements do not access the boundary between 
the lattice and liquid phase of Josephson-vortex system. 
In the future, further improvements of measurements 
would make it possible to investigate the phase diagram 
of Josephson-vortex system including a liquid state by 
using this phenomenon. 

The authors thank X. Hu and M. Machida for useful 
discussions, N. Ishikawa for technical supports of FIB. 
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FIG. 1. An FIB image of a typical sample. Schematic 
drawing of the junction (hatched part) is shown in the inset. 
Magnetic field is applied to the horizontal direction. Defini- 
tion of the junction dimensions is also indicated. 
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FIG. 2. Flow resistance of Josephson vortices as a func- 
tion of magnetic field measured on the sample A (a), B (b), 
and C (c) with the applied currents of 1, 1, and 10 fiA, respec- 
tively. The enlarged figures of the flow-resistance oscillations 
of the sample B and C are shown in the insets of (b) and (c), 
respectively. 
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FIG. 3. I-V characteristics of the sample B at various 
fields from 14 k to 16 kOe at every 20 Oe in 65 K. Each curve 
is shifted by lmV step. Bold lines show I-V curves where the 
resistance is minimum. 
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FIG. 4. Power spectral densities of the oscillations for 
sample A, B and C. Vertical axes are shifted properly. In the 
inset, solid circles show the inverse of H p , which are estimated 
as the positions of the peaks in spectra, as a function of width 
(w). The broken and solid lines represent 1/Hq(w) and its 
twice, respectively. 
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(a) H= n (H /2) (b) H= (n+1) (H /2) 



FIG. 5. Schematic pictures of the configuration of Joseph- 
son- vortex lattice when the magnetic fields are reached to (a) 
n(Ho/2) and (b) (n + l)(i?o/2), where n is integer. 
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